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The  U.S.  Naval  Research  Laboratory  (NRL)  is  developing  nowcast/forecast  software  systems  designed  to 
combine  satellite  ocean  color  data  streams  with  physical  circulation  models  in  order  to  produce  prognos¬ 
tic  fields  of  ocean  surface  materials.  The  Deepwater  Horizon  oil  spill  in  the  Gulf  of  Mexico  provided  a  test 
case  for  the  Bio-Optical  Forecasting  (BioCast)  system  to  rapidly  combine  the  latest  satellite  imagery  of  the 
oil  slick  distribution  with  surface  circulation  fields  in  order  to  produce  oil  slick  transport  scenarios  and 
forecasts.  In  one  such  sequence  of  experiments,  MODIS  satellite  true  color  images  were  combined  with 
high-resolution  ocean  circulation  forecasts  from  the  Coupled  Ocean-Atmosphere  Mesoscale  Prediction 
System  (COAMPS®)  to  produce  96-h  oil  transport  simulations.  These  oil  forecasts  predicted  a  major  oil 
slick  landfall  at  Grand  Isle,  Louisiana,  USA  that  was  subsequently  observed.  A  key  driver  of  the  landfall 
scenario  was  the  development  of  a  coastal  buoyancy  current  associated  with  Mississippi  River  Delta 
freshwater  outflow.  In  another  series  of  experiments,  longer-term  regional  circulation  model  results  were 
combined  with  oil  slick  source/sink  scenarios  to  simulate  the  observed  containment  of  surface  oil  within 
the  Gulf  of  Mexico.  Both  sets  of  experiments  underscore  the  importance  of  identifying  and  simulating 
potential  hydrodynamic  conduits  of  surface  oil  transport.  The  addition  of  explicit  sources  and  sinks  of 
surface  oil  concentrations  provides  a  framework  for  increasingly  complex  oil  spill  modeling  efforts  that 
extend  beyond  horizontal  trajectory  analysis. 
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1.  Introduction 

On  20  April  2010  the  deep-sea  drilling  unit  Deepwater  Horizon 
(DWH)  exploded  leading  to  an  unprecedented  discharge  of  oil  and 
gas  from  the  Macondo  prospect  (~77  km  southeast  of  the  Missis¬ 
sippi  River  Delta;  Fig.  1)  into  the  Gulf  of  Mexico  for  the  following 
86  days.  Estimates  for  the  total  amount  of  oil  released  during  that 
period  range  from  approximately  4.8-8.3  x  10s  L  (Crone  and  Tol¬ 
stoy,  2010;  Leifer,  2010).  This  constitutes  the  largest  accidental 
marine  oil  spill  in  U.S.  waters  (Levy  and  Gopalakrishnan,  2010). 
The  total  economic  impact  of  the  DWH  oil  spill  is  estimated  to 
be  greater  than  US  $8.7  billion  (Sumaila  et  al.,  2012). 

The  unprecedented  scope  of  the  oil  spill  became  obvious  as  sa¬ 
tellite  images  of  the  surface  oil  emerged  in  the  weeks  immediately 
following  the  DWH  blowout.  For  example,  NASA’s  Moderate  Reso¬ 
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lution  Imaging  Spectroradiometer  (MODIS)  sensors  aboard  the 
sun-synchronous  Terra  and  Aqua  satellites  provided  true  color 
images  that  revealed  an  apparent  oil  slick  distribution  contaminat¬ 
ing  >20,000  km2  of  ocean  surface  over  the  course  of  the  oil  spill 
event  (Hu  et  al.,  201 1 ).  Detection  of  the  oil  slick  extent  from  pas¬ 
sive  visible  remote  sensing  is  based  on  the  oil  slick’s  modification 
of  sun  glint  reflectance  (Hu  et  al.,  2009).  Such  detection  methods 
do  not  provide  a  direct  quantitative  assessment  of  oil  concentra¬ 
tion  or  surface  oil  slick  thickness.  Nonetheless,  the  subsequent 
MODIS  and  other  remote  sensing  images  indicated  the  oil  spill 
was  unfolding  as  a  mesoscale  phenomenon  (on  the  spatial  order 
of  ~10-1000  km,  and  temporal  duration  of  weeks  to  months). 

Of  paramount  concern  to  government  agencies,  resource  man¬ 
agers,  and  emergency  responders  during  the  DWH  oil  spill  time 
period  (20  April-15  July  2010)  and  thereafter  was  the  ultimate  fate 
and  potential  landfall  of  the  extensive  offshore  aggregation  of  sur¬ 
face  oil.  Anticipation  of  landfall  required  mobilization  of  extensive 
resources  to  deploy,  for  example,  prophylactic  oil  boom-type  bar- 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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Fig.  1.  Map  of  the  Gulf  of  Mexico  with  Sea  Surface  Height  (SSH)  contours  (positive  solid,  negative  dashed)  provided  by  the  IASNFS  (29  May  2010).  The  Loop  Current  (LC)  and 
an  associated  Loop  Current  Eddy  (LCE)  are  indicated,  inset  and  (b):  detailed  map  of  the  Louisiana  coastal  region  near  the  DWH  site.  Bathymetry  is  indicated  with  dashed  lines. 


riers,  and  to  stage  secondary  defense  supply  stations  in  support  of 
cleanup  efforts  (State  of  Louisiana,  2010). 

The  National  Oceanic  and  Atmospheric  Administration  (NOAA) 
was  the  lead  U.S.  Government  agency  for  oil  slick  trajectory  fore¬ 
casting.  NOAA  provided  nowcasts  of  the  oil  slick  distribution  by 
combining  aircraft  overflights,  satellite  information,  and  in  situ 
observations  (NOAA  OR&R,  2013a).  Forecasts  of  the  oil  slick  distri¬ 
bution  (24,  48,  and  72  h)  were  provided  from  22  April  to  23  August 
2010  (ibid.).  The  forecasting  was  accomplished  via  the  General 
NOAA  Oil  Modeling  Environment  (GNOME)  oil  spill  trajectory 
model  (Zelenke  et  al.,  2012).  The  GNOME  system  ingests  surface 
current  information  from  data  sources  and/or  numerical  ocean  cir¬ 
culation  models  as  well  as  an  initial  oil  contaminant  distribution  to 
project  the  movement  of  these  contaminants.  The  primary  trans¬ 
port  calculation  is  Lagrangian:  i.e.,  surface  oil  is  represented  as  vir¬ 


tual  “particles”  that  are  tracked  over  the  timescale  of  the 
simulation  using  two-dimensional  surface  displacement  calcula¬ 
tions.  This  is  the  common  method  used  in  oil  spill  modeling  (e.g., 
Sotillo  et  al.,  2008)  and  similar  Lagrangian  particle-based  forecast 
methods  were  simultaneously  employed  by  the  research  and  aca¬ 
demic  communities  during  the  DWH  oil  spill  (Liu  et  al.,  2011a; 
Mariano  et  al.,  2011). 

Here  we  present  an  alternative  method  to  two-dimensional 
Lagrangian  oil  trajectory  forecasts.  The  BioCast  system  resolves  a 
fully  three-dimensional  Eulerian  transport  calculation.  These  cal¬ 
culations  do  not  require  presumptions  about  virtual  particles  and 
instead  treat  oil  as  an  idealized  passive  tracer.  Both  types  of  oil  spill 
modeling  approaches  must  make  assumptions  about  the  nature  of 
oil  in  water  that  are  imperfect:  oil  may  behave  as  both  particle 
aggregations  and  dissolved  tracers  depending  on  the  state  of 
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weathering,  dissolution,  and  other  specific  material  properties  of 
the  hydrocarbons  under  consideration  (ASCE,  1996;  Leifer  et  al„ 
2006).  In  a  fundamental  sense,  the  Lagrangian  model  particle  (or 
element)  is  simply  a  point  in  two-dimensional  space  and  its  formal 
representation  of  mass  is  arbitrary.  This  is  in  contrast  to  the  Eule- 
rian  methods  explored  herein:  the  mass  of  oil  in  each  spatial  dis¬ 
cretization  (grid  cell)  defines  the  model’s  state  variable.  This 
allows  for  explicit  calculation  of  three-dimensional  material  trans¬ 
port,  weathering  (transformation)  of  materials,  and  potential  changes 
in  material  buoyancy.  Precedent  for  Eulerian  approaches  to  oil  spill 
modeling  may  be  found  in  Tulloch  et  al.  (201 1 )  and  Maltrud  et  al. 
(2010).  Note  that  in  these  studies  the  tracer  is  described  as  a  generic 
dye,  whereas  herein  we  attempt  to  move  forward  with  an  explicit  oil 
concentration  model.  With  this  increase  in  complexity,  however,  is 
the  associated  disadvantage:  the  modeler’s  dilemma,  i.e„  the  need 
to  parameterize  and  mathematically  represent  processes  that  may 
not  be  not  well  constrained  with  observations  or  experiment. 

Cognizant  of  these  and  other  inherent  complexities,  we  nonethe¬ 
less  seek  to  address  the  remaining  core  problems  posed  by  opera¬ 
tional  oil  transport  forecasting  as  an  oil  spill  response  tool.  First, 
the  methods  required  to  rapidly  combine  satellite-based  estimates 
of  the  oil  spill  distribution  with  state-of-the-art  ocean  circulation 
models  to  produce  an  oil  spill  distribution  forecast  are  evaluated. 
Second,  we  examine  how  the  inherent  reactivity  of  the  contami¬ 
nants  may  impact  the  simulated  distribution  over  time  and  in  con¬ 
trast  to  a  scenario  wherein  only  the  physical  transport  is  considered. 

In  this  paper,  two  Eulerian  transport  approaches  to  oil  spill  sim¬ 
ulation  are  examined  via  numerical  experiment  to  address  these 
aforementioned  problems.  In  the  first  approach  (Section  3),  the 
oil  transport-forecasting  problem  is  examined  in  terms  of  a  passive 
tracer  transported  at  the  ocean’s  surface.  Emphasis  is  placed  on  the 
initial  spatial  distribution  estimated  from  satellite  data  and  the 
evolution  of  this  distribution  over  the  ensuing  96  h.  The  results 
are  examined  with  respect  to  subsequent  observation  of  oil  distri¬ 
bution  and  landfall.  In  the  second  approach  (Section  4),  more  com¬ 
plex  computations  of  potential  hydrocarbon  sources  and  sinks  are 
considered.  These  computations  are  performed  in  the  context  of  a 
longer-term  simulation  (~18days)  of  the  oil  spill  to  address  the 
timescale  of  decay  rate  processes.  Accordingly,  a  regional  ocean 
circulation  model  is  used  to  provide  coastal  surface  current  infor¬ 
mation  as  well  as  simulate  the  interaction  of  the  oil  slick  with 
the  mesoscale  circulation  in  the  open  Gulf  of  Mexico. 

2.  Methods 

BioCast  is  computational  software  that  provides  for  a  rapid 
combination  of  ocean  circulation  model  results  with  a  three- 
dimensional  tracer  transport-reaction  simulation.  The  flow  of 
information  is  thus  very  similar  to  NOAA’s  GNOME  oil  trajectory 
forecasts:  information  on  surface  currents  must  be  combined  with 
some  initial  state  of  the  material  distribution.  BioCast  was  devel¬ 
oped  for  short-term  forecasting  (~24  h)  of  ocean  surface  bio-opti¬ 
cal  properties  as  detected  and  estimated  by  passive  remote  sensing 
methods.  However,  the  software  may  be  applied  to  any  material 
distribution  (such  as  oil)  if  the  initial  state  is  provided. 

The  BioCast  transport  calculation  maps  the  velocity  information 
to  a  three-dimensional  stencil  and  makes  minor  adjustments  to 
constrain  the  velocity  field  to  continuity: 

dw  _du  dv 

~l)z^dx  +  dy  (  ’ 

where  w  is  the  vertical  velocity  and  u,  v  are  the  horizontal  velocities 
in  a  Cartesian  coordinate  frame.  Following  these  adjustments, 
transport  is  calculated  using  first-order  upstream  differences  for 
the  advection  equation  (e.g.,  Smolarkiewicz,  1984)  on  the  three¬ 


dimensional  grid.  First-order  numerical  advection  schemes  are 
inherently  diffusive  (ibid.).  An  analysis  of  the  numerical  diffusion 
inherent  to  our  scheme  yields  horizontal  diffusivity  values  in  the 
range  of  ~50-700  m2  s-1.  Studies  of  natural  horizontal  diffusion 
tend  to  scale  with  the  length  scale  of  the  observations  (Obuko, 
1971).  For  the  length  scales  commensurate  with  the  distribution 
of  oil  slicks  within  our  model  domains  (~50-500  km),  estimates 
of  horizontal  diffusion  are  in  a  similar  range  (~50-1000  m2  s_1; 
Obuko,  1971;  Ledwell  et  al„  1998). 

Thus  the  BioCast  transport  calculation  represents  the  advection- 
diffusion  portion  of  the  advection-diffusion-reaction  problem.  The 
reaction  portion  may  be  modified  in  the  BioCast  software  or  elimi¬ 
nated  entirely  based  on  the  requirements  of  the  problem  and  the  de¬ 
signs  of  the  investigator.  The  reaction  calculations  can  range  from 
simple  decay  rate  constants  to  complex  biogeochemical  models. 
Here  the  reaction  portion  was  modified  to  describe  positive  buoy¬ 
ancy,  and  then  subsequently  modified  to  provide  an  oil  source  term 
and  to  simulate  oil  weathering,  as  explained  in  Section  4. 

In  both  series  of  experiments,  the  initial  state  was  based  on  the 
MOD1S  true  color  imaging  of  the  oil  slick  distribution  on  1 1  May 
2010  (Fig.  2a).  As  mentioned  above,  the  apparent  ocean  surface  dis¬ 
coloration  is  based  largely  on  changes  in  sun  glint  reflectance  due 
to  the  oil  slick’s  presence;  the  varying  angular  dependence  of  sun 
glint  makes  quantitative  retrieval  of  oil  slick  thickness  or  quantity 
from  these  images  very  difficult  (Hu  et  al.,  2009).  Accordingly,  the 
image  was  decomposed  to  develop  a  novel  algorithm  for  determin¬ 
ing  the  spatial  extent  of  the  oil  slick.  Image  pixels  where  oil  is  pre¬ 
sumed  to  be  present  (based  strictly  on  the  apparent  contrast  with 
the  surrounding  open  ocean  image  pixels)  are  analyzed  for  the 
scaled  red  (R)  and  green  (G)  image  values  (ranging  from  1  to 
255).  This  is  repeated  (~10  times)  to  develop  an  approximate 
range  of  values  for  apparent  oil-influenced  surface  water  discolor¬ 
ation  in  the  true  color  image.  Once  a  set  of  thresholds  is  estab¬ 
lished,  all  of  the  oil-containing  pixels  are  identified  via 
automated  image  processing  software.  This  procedure  must  be  re¬ 
peated  and  adjusted  for  any  new  RGB  image  because  the  RGB  true 
color  data  processing  will  render  different  scaled  RGB  values  based 
on  the  amount  of  sun  glint  reflectance  present  in  the  raw  satellite 
data.  A  similar  procedure  was  also  used  to  remove  the  presence  of 
clouds,  again  based  on  the  RGB  values  where  clouds  were  pre¬ 
sumed  to  be  present.  Image  pixels  identified  as  oil  (Fig.  2b)  were 
then  mapped  to  the  corresponding  latitude  and  longitude  coordi¬ 
nates  of  the  model  domain,  and  thereby  provided  a  starting  place 
for  forward  integration  of  the  transport  computations  (Fig.  3). 
Additional  steps  to  initialize  the  surface  oil  concentration  on  a 
more  quantitative  basis  are  explained  in  Section  4. 

There  is  uncertainty  in  this  initial  surface  oil  distribution.  Other 
concurrent  satellite  analyses  depict  additional  surface  oil  west  and 
north  of  our  estimate  (Hu,  2010),  and  potentially  smaller  oil  slicks 
detached  from  the  main  bolus  near  the  DWH  site  (NOAA/NESDIS, 
2013).  We  note  that  synthetic  aperture  radar  (SAR)  sensors  also  pro¬ 
vided  satellite-based  estimates  of  surface  oil  slick  locations  during 
the  event  that  may  have  been  substantially  different  from  MODIS 
sun  glint-based  analyses  (Walker  et  al.,  2012).  Future  work  will 
aim  towards  a  more  comprehensive  assimilation  of  satellite  infor¬ 
mation  with  associated  uncertainty  estimates  into  oil  spill  models. 

An  initial  set  of  numerical  experiments  was  performed  using 
the  Coupled  Ocean-Atmosphere  Mesoscale  Prediction  System 
(COAMPS®),  a  nested  modeling  system  developed  at  the  Naval  Re¬ 
search  Laboratory  that  allows  for  a  two-way  exchange  of  informa¬ 
tion  between  the  atmospheric  and  oceanic  forecasting 
components.  The  nonhydrostatic  atmospheric  COAMPS  model 
component  (Hodur,  1997)  is  the  operational  mesoscale  forecasting 
system  for  the  U.S.  Navy.  The  hydrostatic  Navy  Coastal  Ocean  Mod¬ 
el  (NCOM;  Martin,  2000;  Barron  et  al.,  2004;  Kara  et  al.,  2006) 
served  as  the  oceanic  component.  NCOM  is  the  main  regional 
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(a) 


Fig.  2.  MODIS  data  acquired  11  May  2010, 18:55  UTC  (250  m  horizontal  resolution)  and  processed  as  a  true  color  image  is  shown  in  plate  (a).  In  (b)  the  apparent  position  of 
the  surface  oil  slicks  are  extracted  from  the  image,  as  explained  in  the  text. 


oceanographic  forecasting  model  for  the  U.S.  Naval  Oceanographic 
Office  (Rhodes  et  al.,  2002). 

The  atmospheric  and  oceanic  model  coupling  was  designated 
via  the  upper-most  oceanic  model  grid  cell  temperature  and  the 
lowest  grid  cell  atmospheric  model  variables  (temperature,  humid¬ 
ity,  wind  velocity,  pressure,  and  radiative  fluxes).  At  a  6-min  cou¬ 
pling  interval,  bulk  fluxes  of  heat  energy  exchange  were 
calculated  following  the  Coupled  Ocean-Atmosphere  Response, 
version  3  (COARE  3.0)  scheme  (Fairall  et  al.,  1996).  Further  details 
of  the  COAMPS  modeling  components  are  listed  in  Small  et  al. 
(2012).  Verification  and  validation  of  the  COAMPS  forecasting  sys¬ 
tem  may  be  found  elsewhere  (Doyle  et  al.,  2009;  Small  et  al., 
2012);  here  we  focus  on  how  the  forecasts  of  surface  currents  from 
COAMPS  may  be  used  by  the  BioCast  system  to  forecast  surface  oil 
transport.  Only  the  “true”  hourly  surface  current  forecast  fields  for¬ 
ward  from  the  analysis  time  (the  initial  state  on  11  May  2010) 
were  used.  This  means  that  the  surface  current  velocities  were  gen¬ 
uine  forecasts  of  marine  conditions  from  the  modeling  system;  i.e., 
no  data  assimilation  of  atmospheric  or  oceanic  data  beyond  the 
analysis  time  occurred. 

The  estimated  surface  oil  distribution  from  the  11  May  2010 
MODIS  image  was  used  to  initialize  a  relative  oil  concentration 
state  variable.  The  oil  concentration  was  treated  as  a  passive  tracer 
(physical  transport/no  biological-chemical  reactions)  with  no 


additional  sources  beyond  the  initialization  field.  Hence  the  con¬ 
servation  equation  may  be  simply  expressed  as: 


dRC 

Ijf 


/  du  dv 
V9x  +  dy  + 


dw 

~dz 


RC  +  B 


(2) 


The  state  variable,  RC,  is  a  relative  concentration  of  oil.  This  va¬ 
lue  was  initialized  as  100  where  the  MODIS  true  color  threshold- 
based  algorithm  suggested  the  presence  of  oil. 

The  transport  calculation  treats  the  surface  oil  as  a  dissolved 
tracer.  This  permits  downwelling  (downward  vertical  advection) 
as  well  as  diffusion  into  grid  cells  beneath  the  surface.  Whereas 
this  may  indeed  be  the  fate  of  some  dissolving  or  emulsified  hyrd- 
ocrabons,  a  positive  buoyancy  term  (B)  was  nonetheless  added  to 
the  calculation  (Eq.  (2))  to  force  the  simulated  hydrocarbons  back 
into  the  surface  grid  cell.  This  calculation  does  not  permit  down¬ 
ward  vertical  penetration  of  simulated  oil  but  it  will  still  allow 
for  surface  convergence  or  divergence  (dispersion)  of  materials. 

Once  again,  this  is  in  contrast  to  Langrangian  methods.  Physical 
dispersion  cannot  be  explicitly  defined  for  a  single  point  in  space 
subject  to  a  horizontal  displacement  calculation.  Statistical  meth¬ 
ods  must  be  employed  to  apply  a  probabilistic  modification  of 
the  trajectory;  e.g.,  the  random  walk  method  (Proehl  et  al., 
2005).  Veracity  of  these  statistical  methods  generally  improves 
with  an  increase  in  the  number  of  representative  Lagrangian 
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Fig.  3.  The  initial  estimate  of  surface  oil  slick  distribution  extracted  from  the  MOD1S  image  is  mapped  to  the  inner  ocean  model  domain  of  COAMPS  (500  m  horizontal 
resolution)  for  the  beginning  of  an  oil  trajectory  forecast.  The  oil  forecast  output  increment  (+1  h)  is  shown  for  1 1  May  2010,  2000  UTC.  RC  is  the  Relative  Concentration, 
scaled  to  the  initial  oil  distribution  estimate,  RC  =  100.  The  dashed  line  indicates  RC=  1,  the  solid  line  begins  the  contours  at  RC  =  5. 


particles  or  alternatively,  an  ensemble  of  trajectory  model  simula¬ 
tions  (Brickman  and  Smith,  2002).  As  the  number  of  trajectory  par¬ 
ticles  tracked  (or  ensembles)  increases,  however,  so  does  the 
computational  expense.  Finally,  one  arrives  near  the  impractical 
end  of  that  continuum  and  may  elect  to  instead  perform  a  Eulerian 
computation  that  explicitly  calculates  the  mass  flux  of  distributed 
materials  in  a  single  iteration.  The  disadvantage  now  is  that  the 
representation  of  the  oil’s  mass  as  uniformly  distributed  over  the 
discrete  spatial  resolution  of  the  model  may  tend  to  result  in  overly 
dispersive  transport.  Here,  the  vertical  dispersion  is  eliminated  via 
a  buoyancy  restoring  term.  Horizontal  dispersion  remains.  A  con¬ 
sequent  criticism  of  this  Eulerian  framework  specific  to  oil  spill 
modeling  is  that  “it  is  practically  impossible  to  detect  exactly  the 
oil  spill  boundaries  in  a  specific  moment”  (Lonin,  1999).  This  is 
in  contrast  to  a  spatial  distribution  of  Lagrangian  elements  that 
provides  a  very  discrete  oil  spill  boundary.  As  a  practical  matter, 
however,  this  criticism  is  easily  addressed.  One  approach  is  to  scale 
the  Eulerian  tracer  field  to  the  initial  source  concentration  (as  in 
Maltrud  et  al.,  2010;  and  here,  Section  3).  Another  approach  is  to 
define  the  horizontal  oil  spill  boundary  using  a  lower  limit  of 
detection,  or  threshold  value  (as  in  Section  4). 

3.  COAMPS-based  forecast  results 

The  24-h  forecast  shows  the  lateral  spreading  of  the  initial  rel¬ 
ative  concentration  (RC)  field.  The  northwest  corner  of  the  oil  slick 
is  initiating  contact  with  the  Mississippi  River  Delta  (MRD;  Fig.  4a). 
The  48-h  forecast  indicates  this  oil  is  being  transported  clockwise 
around  the  Southwest  Pass  of  the  MRD  and  initiating  landfall  on 
the  southern  coast  of  Plaquemines  Parish,  Louisiana  and  towards 
Barataria  Bay  (Fig.  4b).  Some  —74  h  forward  into  the  forecast  per¬ 
iod,  this  clockwise  conduit  around  the  MRD  funnels  increasing 


amounts  of  the  initial  surface  oil  distribution  into  the  Louisiana 
Bight  to  make  significant  landfall  along  the  outer  islands  of  Bara¬ 
taria  Bay,  including  Grand  Isle  and  southwest  towards  Port  Four- 
chon  (Fig.  5a).  At  the  conclusion  of  the  forecast  period  (96  h;  the 
full  sequence  is  provided  in  Animation  1),  this  pattern  persists 
and  much  of  the  remaining  oil  from  the  initial  distribution  is  being 
transported  westward  into  an  apparent  onshore/offshore  bifurca¬ 
tion  in  the  surface  oil  distribution  (Figs.  5b).  Smaller  amounts  of 
oil  have  also  been  transported  northwest  towards  the  Chandeleur 
Islands  and  into  Breton  Sound.  Comparatively,  however,  a  far  larger 
amount  of  the  initial  oil  is  transported  into  the  Louisiana  Bight  to 
ultimately  make  landfall  at  coastal  Louisiana  locations  west  of 
the  MRD. 

The  oil  transport  patterns  are  explained  by  the  concurrent  sur¬ 
face  current  forecasts  obtained  from  COAMPS  (Fig.  6a).  Large  veloc¬ 
ity  surface  currents  (—1.3  m  s_1)  are  moving  clockwise  around  the 
MRD.  This  circulation  feature  is  coherent  and  well-established 
by  approximately  48  h  and  persists  through  the  remainder  of  the 
forecast  period  (Fig.  6b).  There  is  also  a  bifurcation  in  the  surface 
flow  —45  km  south  of  the  MRD  that  explains  the  apparent  off¬ 
shore/onshore  divergence  in  the  simulated  surface  oil  patterns 
(Fig.  6a  and  b). 

This  forecast  of  oil  transport  was  qualitatively  accurate.  Oil  from 
the  DWH  spill  was  observed  making  landfall  in  the  vicinity  of  Port 
Fourchon  on  14  May  2010  (Schmidt,  2010).  Ground  observations 
reported  by  the  Shoreline  Cleanup  Assessment  Technique  (SCAT) 
teams  indicate  initial  landfall  of  oil  on  14  May  2010  along  the  Lou¬ 
isiana  coast  from  Port  Fourchon  to  Grand  Isle  (NOAA  ORSiR,  2013b). 
Heavier  amounts  of  surface  oil  were  sighted  in  the  vicinity  of 
Grand  Isle,  Louisiana  on  20  May  2010  (Rioux,  2010).  By  23  May 
2010,  SCAT  data  indicate  heavy  landfall  of  oil  occurring  from  Port 
Fourchon  to  Barataria  Bay.  This  was  followed  by  some  of  the  most 
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Fig.  4.  (a)  Oil  transport  forecast  for  12  May  2010, 1900  UTC,  and  (b)  oil  transport  forecast  for  13  May  2010, 1900  UTC.  The  dashed  line  indicates  RC=  1,  the  solid  line  begins 
the  contours  at  RC  =  5. 


significant  landfall  of  surface  oil  associated  with  the  DWH  event 
(OSAT-2,  201 1 ).  The  salt  marshes  of  Barataria  Bay  were  also  some 
of  the  most  severely  oiled  coastal  habitats  (Michel  et  al.,  2013;  Zen- 
gel  and  Michel,  2013). 


In  addition  to  the  ground  observations,  the  forecast  results  are 
confirmed  by  concurrent  analysis  of  satellite  imagery.  NOAA 
National  Environmental  Satellite,  Data  and  Information  Service 
(NESDIS)  satellite  composite  analysis,  which  incorporates  MODIS 
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Fig.  5.  (a)  Oil  transport  forecast  for  14  May  2010,  2100  UTC,  and  (b)  oil  transport  forecast  for  15  May  2010, 1900  UTC.  The  dashed  line  indicates  RC  =  1,  the  solid  line  begins 
the  contours  at  RC  =  5. 


data,  SAR  data,  and  other  sensors  (NOAA/NESD1S,  2013),  verifies 
the  movement  of  large  oil  slicks  into  the  Louisiana  Bight  on  17 
May  2010  (Fig.  7a  and  b).  The  20  May  2010  analysis  suggests  the 
conduit  around  the  Southwest  Pass  was  indeed  persistent.  Addi¬ 
tional  oil  more  than  45  km  directly  south  of  the  MRD  would  also 


support  the  bifurcation  in  surface  currents  depicted  in  the  COAMPS 
forecast  and  manifest  in  the  simulated  oil  distributions.  The  23 
May  2010  NOAA/NESD1S  analysis  indicates  surface  oil  in  Breton 
Sound  and  around  the  Chandeleur  Islands  (Fig.  7c).  SCAT  data 
confirm  concurrent  landfall  in  Chandeleur  Sound.  The  23  May 
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Fig.  6.  (a)  COAMPS  surface  current  velocity  field  forecast  for  14  May  2010,  1700  UTC,  and  (b)  COAMPS  surface  current  velocity  field  forecast  for  14  May  2010, 1700  UTC. 


NOAA/NESDIS  analysis  also  depicts  oil  entering  Barataria  and  Ter¬ 
rebonne  Bays  (Fig.  7c). 

Given  the  qualitative  agreement  between  forecast  and  observa¬ 
tions,  it  is  probable  that  the  forecast  surface  current  fields  have 
some  fidelity  to  genuine  surface  currents  between  11  and  19 
May  2010.  However,  our  forecasts  based  on  the  11  May  2010  ini¬ 
tialization  apparently  accelerated  the  landfall  of  significant  surface 
oil  slicks  upon  Grand  Isle,  Louisiana  and  vicinity  to  14  May  2010, 
whereas  observations  suggest  heavy  landfall  of  oil  did  not  truly 
commence  until  ~19-20  May  and  thereafter.  The  SCAT  data  record 
of  landfall  in  these  areas  on  14  May  2010  is  documented  as  “very 
light,”  i.e.,  consisting  of  isolated  pockets  of  tar  balls  and  scattered 


emulsified  oil  aggregations.  Other  ground  observations  verify  this 
description  (Schmidt,  2010).  More  severe  categories  of  land  surface 
oiling  appear  to  commence  in  the  SCAT  record  around  20  May 
2010.  Part  of  the  temporal  discrepancy  between  our  forecast  and 
observations  may  be  due  to  oil  weathering  and  the  application  of  dis¬ 
persants.  These  processes  were  not  represented  in  these  COAMPS- 
based  forecasting  experiments.  Another  possibility  is  that  landfall  of 
oil  is  a  process  that  is  distinct  from  shoreward  progression  and  its 
simulation  requires  model  parameterizations  of  winds,  surface  waves, 
and  littoral  tidal  processes  below  the  spatial  resolution  of  our  models. 

Another  potential  source  of  the  temporal  mismatch  may  be 
associated  with  the  development  and  intensification  of  a  buoyancy 
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Mississippi  River  discharge  data  (Tarbert  Landing,  MS)  provided 
by  the  United  States  Army  Corps  of  Engineers  indicate  below  his¬ 
torical  values  from  20  April  to  11  May  2010  (Fig.  8).  This  may 
partly  explain  the  temporal  mismatch  between  simulated  and  ob¬ 
served  oil  slick  landfall:  the  simulated  buoyancy  current  was  well 
established  by  13  May  whereas  the  true  currents  may  have  been 
less  intense  until  sustained  freshwater  discharge  out  of  the  MRD 
was  sufficient  to  accumulate  a  substantially  larger  buoyancy 
plume. 

This  buoyancy  current  is  a  recurrent  and  characteristic  feature 
in  this  area  (Rouse  et  al„  2005).  It  is  thus  highly  probable  that 
any  oil  spill  in  the  vicinity  of  the  MRD  will  make  landfall  along 
Grand  Isle,  Louisiana  and  the  adjacent  coastal  sections  of  the  Lou¬ 
isiana  Bight.  This  landfall  would  encompass  Plaquemines,  lower 
Jefferson,  and  LaFourche  Parishes  and  would  potentially  propagate 
farther  west  to  Terrebonne  Bay  (see  Fig.  1).  Emergency  managers 
and  government  agencies  should  be  cognizant  of  this  probability. 


Fig.  7.  Composite  satellite  analysis  of  potential  surface  oil  location  obtained  from 
the  NOAA/NESD1S  archive  for  (a)  17  May  2010,  (b)  20  May  2010,  and  (c)  23  May 
2010. 


current  along  the  MRD  and  the  upper  Louisiana  Bight.  Generally, 
buoyant  spreading  of  low  salinity  water  from  a  river  mouth/delta 
or  estuary  in  the  Northern  Hemisphere  will  propagate  with  land 
to  the  right  (looking  down  current)  (Simpson,  1997).  Along  the 
Louisiana  Bight  and  Louisiana-Texas  coasts,  this  recurrent  coastal 
circulation  feature,  augmented  by  southeasterly  winds,  is  known 
as  the  Louisiana  Coastal  Current  (LCC)  (Wiseman  and  Dinnel, 
1988). 


4.  Regional  source/sink  experiments 

4.1.  Velocity  fields  and  oil  initialization 

The  first  series  of  experiments  represented  the  surface  oil  as  a 
buoyant  tracer  and  focused  on  a  96-h  forecast  within  the  inner  nest 
(500  m  horizontal  resolution)  of  a  nested  ocean  modeling  domain. 
For  a  second  set  of  numerical  experiments,  the  domain  was  ex¬ 
panded  to  include  the  entire  Gulf  of  Mexico  and  incorporate  results 
from  a  regional  ocean  circulation  model.  The  Intra-Americas  Seas 
Nowcast/Forecast  System  (1ASNFS;  Ko  et  al.,  2003)  provided  regio¬ 
nal  (~3  km  horizontal  resolution)  flow  fields  for  integration  into 
the  BioCast  system.  IASNFS  is  a  regional  application  of  NCOM. 
The  Navy  Operational  Global  Atmospheric  Prediction  System  (NO¬ 
GAPS)  provided  atmospheric  surface  forcing  (Rosmond,  1992). 

These  series  of  experiments  are  not  true  forecasts  in  the  sense 
that  the  ocean  circulation  model  results  are  taken  from  the  analysis 
fields.  The  term  “analysis  fields”  refers  to  the  assimilation  of  satel¬ 
lite  data  into  the  modeling  system  via  the  Modular  Ocean  Data 
Assimilation  System  (MODAS)  (Fox  et  al.,  2002).  MODAS  assimi¬ 
lates  remotely-sensed  sea  surface  temperature  (SST)  and  sea  sur¬ 
face  height  (SSH)  data  that  have  been  optimally  interpolated 
(Bretherton  et  al.,  1976)  onto  a  two-dimensional  grid.  Potential 
subsurface  temperature  departure  from  a  long-term  climatology 
(U.S.  Navy  Master  Ocean  Observation  Database  -  MOODS)  is  then 
calculated  based  on  regression  coefficients  that  derive  subsurface 
temperature  from  SSH  and  SST.  The  result  is  a  synthetic  three- 
dimensional  temperature  field.  The  combined  SST  and  SSH 
assimilation  provides  fidelity  to  the  mesoscale  dynamics  in  the 
Gulf  of  Mexico,  which  is  critical  to  forecasting  the  regional-scale 


Fig.  8.  Discharge  data  (solid  line)  for  Mississippi  River  at  Tarbert  Landing  (United 
States  Army  Corps  of  Engineers)  are  shown  in  cubic  meters  per  second.  The 
asterisks  (*)  indicate  the  climatological  monthly  means  used  for  the  COAMPS 
simulations. 
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Fig.  9.  (a)  Regional  oil  model  result  for  17  May  2010, 1800  UTC.  SSH  contours  (positive  solid,  negative  dashed)  provided  by  the  corresponding  IASNFS  fields,  (b)  MOD1S  true 
color  image  acquired  on  17  May  2010,  1640  UTC.  (c)  The  estimate  of  visible  oil  is  mapped  to  the  model  domain  and  shown  in  black. 


circulation.  MODAS  synthetics  have  been  previously  used  to  exam¬ 
ine  biophysical  patterns  in  the  Gulf  (Jolliff  et  al.,  2008).  The  ar¬ 
chived  IASNFS  analysis  fields  are  more  properly  considered 
“hindcasts.” 

Oil  initialization  was  again  based  on  the  MODIS  11  May  2010 
image.  An  accurate  quantitative  estimate  of  surface  oil  concen¬ 
tration  based  solely  on  sun  glint  reflectance  in  a  satellite  image 
is  not  presently  feasible.  However,  it  is  reasonable  to  presume 
there  is  some  minimum  threshold  of  oil  presence  at  the  ocean 
surface  that  must  be  reached  before  any  detection  with  passive 
visible  remote  sensing  techniques  may  occur.  Modification  of 
sun  glint  reflectance  by  surface  oil  suggests  the  presence  of  oil 
in  sufficient  thickness  to  suppress  short  surface  waves  (Adamo 
et  al.,  2009).  With  respect  to  operational  monitoring,  an  oil 
“slick”  is  defined  as  oil  of  sufficient  thickness  to  dampen  surface 
waves  (NOAA  OR&R,  2012).  Based  on  charts  adapted  from  the 
Bonn  Agreement  Oil  Appearance  Codes  (BAOAC)  (ibid.),  minimum 
satellite  detected  oil  slick  thickness  is  estimated  to  be  2.5  pm. 


Note  that  this  is  different  from  the  minimum  thickness  visible 
to  the  human  eye.  Here  the  estimate  is  focused  on  the  minimum 
thickness  for  MODIS  sun  glint-contaminated  images  of  surface  oil. 
Using  a  standard  reference  density  for  Texas  crude  oil  (873  kg  itT3), 
the  model  is  initialized  at  2180  mg  oil  per  irr2  of  ocean  surface  for 
those  grid  cells  where  we  presume  the  presence  of  oil  from  the 
MODIS  true  color  image  (Fig.  2b). 

The  initial  concentration  values  are  determined  by  dividing 
the  initial  per  unit  area  value  (2180  mg  oil  irr2)  by  the  depth 
of  the  surface  grid  cell.  As  before,  a  positive  buoyancy  restoring 
term  maintains  the  oil  in  the  model’s  surface  grid  cells.  The 
model  results  are  converted  back  to  a  per  unit  surface  area  basis 
for  analysis  (Fig.  9a).  In  reality,  some  hydrocarbons  may  be  dis¬ 
solved  whereas  much  remains  at  the  surface  to  form  slicks  and 
sheen.  If  one  assumes  all  of  the  simulated  oil  per  unit  area  in 
the  model  is  at  the  surface  in  the  form  of  a  surface  slick,  then 
the  thickness  of  the  slick  (or  sheen)  may  be  calculated  using 
the  reference  density. 
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Fig.  10.  (a)  Regional  oil  model  result  for  19  May  2010, 1200  UTC.  (b)  Regional  oil  model  result  for  22  May  2010,  0600  UTC. 


4.2.  Source  term 

Here  we  consider  a  source  term  based  on  oil  apparent  at  the 
ocean’s  surface.  Oil  entering  the  water  at  depth  (~1500m)  was 
not  explicitly  modeled.  There  are  likely  many  processes  impacting 
oil  injected  at  depth  that  curtail  its  subsequent  appearance  at  the 
surface  (Socolofsky  et  al„  2011;  Joye  et  al.,  2011).  For  this  reason, 
an  initial  source  term  was  added  at  the  surface  of  the  DWH  site: 
32.2  Ls_1,  or  approximately  17,500  barrels  per  day  (BPD).  This  esti¬ 
mate  was  based  on  the  surface  mass  balance  estimates  provided  by 
the  National  Incident  Command,  Flow  Rate  Technical  Group 
(13,000-22,000  BPD)  (McNutt  et  al.,  2011).  Once  again  applying  a 
standard  reference  density  for  Texas  crude,  a  mass  flux  of 
28  kg  oil  s_1  is  added  at  the  grid  cell  encompassing  the  DWH  site 
within  the  model  domain. 


4.3.  Sink  term 

A  simple  first-order  decay  rate  estimate  provided  a  sink  term  to 
account  for  evaporation,  weathering  and  removal  processes  other 
than  physical  transport.  Evaporation  of  crude  oil  has  been  shown 
to  follow  simple  decay  rate  kinetics  (Fingas,  1995),  and  evapora¬ 
tion  is  an  important  process  with  respect  to  mass  balance  of 
surface  oil  (National  Research  Council,  2003).  Given  that  “Missis¬ 
sippi  Canyon  252  crude  oil”  (Belore  et  al.,  2011)  will  experience 
45%  loss  from  surface  evaporation  after  2  weeks  (ibid.),  and  given 
a  first-order  decay  rate  relation: 
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Fig.  11.  (a)  Regional  oil  model  result  for  29  May  2010,  0600  UTC.  (b)  Regional  oil  model  result  for  the  NL  simulation  29  May  2010,  0600  UTC. 


the  decay  rate  constant  (r)  is  4.9  x  10  7  s  \  The  decay  term  in  the  dC 
model  is  then  simply:  dt 


fdu  dv  dw\ 
[dx  +  dy+^z) 


C+B 


+  a(i,j)  -  rC 


(5) 


C  =  C0e-n  (4) 

The  processes  contributing  to  the  weathering  and  removal  of 
surface  oil  are  multifaceted  and  complex.  Biodegradation  may  re¬ 
move  some  lower  molecular  weight  hydrocarbons  from  the  bulk 
crude  oil  on  much  shorter  timescales,  whereas  higher  molecular 
weight  compounds  may  be  much  more  recalcitrant  to  biodegrada¬ 
tion  and  weathering  processes  (Atlas  and  Hazen,  2011).  The  point 
of  this  numerical  simplification  is  simply  to  establish  the  timescale 
of  the  overall  oil  slick  degradation.  Obviously,  some  components  of 
the  crude  would  require  (r)  values  in  Eq.  (4)  much  larger  or  smaller 
than  4.9  x  10_7s_1;  however,  incorporating  this  higher  level  of 
complexity  into  the  simulation  requires  significant  expansion  of 
the  state  variable  space,  and  hence  demands  for  additional  detailed 
knowledge  of  source  concentrations  and  chemical  composition. 

Given  these  considerations,  the  conservation  equation  for  sur¬ 
face  oil  is  given  as 


where  the  change  in  surface  oil  (C;  mg  oil  per  m3  ocean  surface) 
with  respect  to  time  is  the  transport  calculation  plus  ( B ),  the  buoy¬ 
ancy  calculation,  and  the  source/sink  terms.  The  source  term  (a)  is 
zero  everywhere  except  the  surface  location  (i,  j)  of  the  DWH  site, 
and  r  is  the  universal  decay  rate  constant. 

4.4.  Regional  source/sink  simulation  results 

Large  portions  of  the  DWH  oil  slicks  are  simulated  to  entrain 
into  the  outer  edge  of  the  Loop  Current  (Animation  2),  as  indicated 
by  the  IASNFS  model’s  sea  surface  height  contours  (Fig.  9a).  This 
large  anti-cyclonic  feature  in  the  northern  Gulf  is  almost  pinching 
off  from  the  Loop  Current  to  form  a  warm-core  eddy.  The  large 
extension  of  oil  slick  into  the  Gulf  simulated  on  17  May  is  qualita¬ 
tively  similar  to  the  MOD1S  true  color  image  captured  on  1 7  May 
(Fig.  9b  and  c).  It  is  reasonable  to  conclude  that  such  oil  features 


96 


J.K.  Jolliffet  al./Ocean  Modelling  75  (2014)  84-99 


Fig.  12.  (a)  Regional  oil  model  result  for  29  May  2010,  0600  UTC  using  the  (4 r)  decay  rate  sensitivity  simulation  is  contoured  as  in  Fig.  11.  The  contour  of  the  NL  simulation 
10  mg  oil  m  2  surface  isopleth  is  shown  in  red.  (b)  Sensitivity  analysis  for  surface  oil  concentrations  indicated  in  (a).  The  surface  concentration  values  for  each  different 
simulation  (varying  values  of  r)  is  normalized  by  the  corresponding  value  in  the  NL  simulation  and  expressed  as  a  percentage. 


(once  transported  into  the  Loop  Current)  will  likely  transit  out  of 
the  Gulf  and  into  the  Florida  Straits.  Indeed,  there  was  speculation 
supported  by  trajectory  model  evidence  during  the  oil  spill  that 
this  may  potentially  occur  (Nelson,  2010).  As  the  oil  spill  pro¬ 
ceeded,  however,  no  significant  surface  oil  transport  out  of  the  Gulf 
of  Mexico  was  observed  (Liu  et  al„  2011b). 

The  hindcast  simulations  suggest  two  main  reasons  for  this  fail¬ 
ure  to  egress  the  Gulf.  First,  much  of  the  simulated  oil  initially 
transported  offshore  appears  to  recirculate  within  a  cyclonic  eddy 
associated  with  the  Loop  Current.  Walker  et  al.  (2012)  refer  to  this 
feature  as  a  Loop  Current  frontal  eddy  and  document  its  evolution 
during  the  oil  spill.  In  our  simulation,  the  remaining  surface  oil 
does  not  genuinely  entrain  into  the  outer  Loop  Current  until  such 
time  as  this  larger  anti-cyclonic  circulation  feature  has  finally  de¬ 
tached  to  form  a  Loop  Current  Eddy  (LCE).  A  secondary  and  aug¬ 
menting  mechanism  of  Gulf  containment  in  our  simulation  is  the 
decay  rate  term  that  significantly  degrades  simulated  surface  oil 
and  thereby  reduces  its  horizontal  extent. 


Regarding  the  LCE,  such  anti-cyclonic  circulation  features  fre¬ 
quently  detach  from  the  Loop  Current  and  propagate  westward 
in  the  Gulf  (Leben  and  Born,  1993),  and  such  events  are  often  asso¬ 
ciated  with  the  appearance  of  cyclonic  circulation  features  (Biggs 
et  al.,  1996).  Two  cyclonic  circulation  features  are  evident  on  17 
May  (Fig.  9a):  one  where  the  inchoate  LCE  appears  to  be  detaching 
from  the  Loop  Current,  and  another  at  the  top  of  the  LCE  where  the 
leading  edge  of  the  simulated  oil  plume  appears  to  be  entering  a 
convergent  circulation  feature.  In  the  simulation,  the  surface  oil 
is  recirculating  within  this  feature  (Fig.  9a)  whereas  in  the  MODIS 
image  the  oil  “trail”  extending  into  the  Gulf  appears  to  be  just 
beginning  a  turn  towards  the  northeast  at  its  apparent  terminus 
(Fig.  9c).  This  feature  is  also  depicted  in  the  17  May  NOAA/NESD1S 
analysis  (Fig.  7a). 

Some  ~2  days  forward  in  the  simulation,  the  offshore  oil  still 
appears  to  be  circulating  in  the  cyclonic  frontal  eddy  northeast  of 
the  main  LCE  (Fig.  10a).  Corroborating  evidence  of  this  surface 
entrapment  of  oil  within  a  convergent  circulation  pattern  is  shown 


J.K.  Jolliff  et  al./Ocean  Modelling  75  (2014)  84-99 


97 


in  the  20  May  satellite  analysis  image  (Fig.  7b).  This  general  off¬ 
shore  pattern  of  surface  oil  persists  into  22  May  with  the  addition 
of  some  trace  amounts  of  oil  penetrating  around  the  periphery  of 
the  LCE,  which  has  now  finally  detached  from  the  Loop  Current 
(Fig.  10b).  Note  that  the  contour  intervals  in  the  oil  plots  terminate 
at  10  mg  oil  m  2  of  ocean  surface.  Given  the  assumptions  pre¬ 
sented  in  Section  4.1,  this  would  correspond  to  surface  oil  sheen 
of  approximately  0.01  pm  in  thickness.  This  is  below  the  minimum 
threshold  of  surface  oil  appearance  in  the  BAOAC  charts  (0.04  pm; 
NOAA  OR&R,  2012).  Thus  only  a  trace  amount  of  oil  appears  to  fi¬ 
nally  transit  around  the  periphery  of  the  LCE. 

To  elucidate  the  potential  role  of  degradation/weathering  in  the 
simulated  oil  distribution,  a  second  experiment  was  performed 
wherein  the  decay  rate  constant  (r)  was  set  to  zero:  a  no  loss 
(NL)  simulation.  All  other  aspects  of  the  simulation  were  identical 
to  the  initial  case.  In  the  final  output  frame  of  both  simulations, 
17.5  days  after  the  initial  start-up,  the  same  overall  spatial  distri¬ 
bution  is  evident  (Fig.  11).  As  earlier,  some  offshore  oil  is  recircu¬ 
lating  in  a  cyclonic  frontal  eddy  northeast  of  the  LCE  and  some 
trace  of  oil  is  beginning  to  circulate  around  the  outer  edge  of  the 
LCE.  The  leading  edge  of  this  oil  plume  extends  approximately 
360  km  farther  in  the  NL  simulation  (Fig.  lib).  Elsewhere,  the  NL 
simulation  depicts  larger  amounts  of  oil  at  the  surface  where  oil 
is  present. 

The  role  of  simulated  decay  in  the  surface  oil  distribution  is  ex¬ 
plored  further  with  a  sensitivity  analysis  of  the  decay  rate  constant. 
The  horizontal  extent  of  surface  oil  in  NL  simulation  (defined  by 
the  10  mg  itT2  surface  oil  isopleth)  is  reduced  by  42%  when  the  de¬ 
cay  rate  constant  (r)  is  present  and  increased  by  a  factor  of  four 
(Fig.  12a).  Surface  oil  concentration  is  also  evaluated  at  three  dif¬ 
ferent  locations  for  the  final  frame  of  the  simulation  (17.5  days): 
(1)26  km  southwest  of  the  DWH  site,  (2)  in  the  center  of  the  cyclo¬ 
nic  frontal  eddy  (277  km  from  DWH),  and  (3)  along  the  outer  edge 
of  the  LCE  (464  km  from  DWH;  Fig.  12).  The  far  field  sites  (2  and  3) 
are  significantly  impacted  by  the  decay  rate  (Fig.  12b).  If  (r)  is  in¬ 
creased  by  a  factor  of  four  then  the  final  concentration  for  both 
sites  is  reduced  to  below  11%  of  the  corresponding  NL  simulation 
value.  These  results  reveal  a  very  large  sensitivity  for  the  far  field 
sites  between  0.25r  and  4r  (~82%  to  <11%).  This  corresponds  to  a 
half-life  decay  of  65.5  down  to  4.1  days.  In  contrast,  the  near  field 
site  (1 )  concentrations  are  all  greater  than  60%  of  the  NL  value  over 
the  entire  range  of  (r)  values. 


5.  Discussion  and  conclusions 

These  regional  Gulf  of  Mexico  oil  spill  simulations  demonstrate 
how  the  southern  Florida  coastline  was  spared  contact  with  any 
significant  bolus  of  surface  oil  due  to  the  fortuitous  arrangement 
of  mesoscale  circulation  features  and  the  subsequent  detachment 
of  a  warm-core  eddy  from  the  Loop  Current.  Had  this  not  occurred, 
however,  our  simulations  suggest  that  the  weathering  and  decay  of 
the  surface  oil  may  have  mitigated  any  potential  impact.  We  note 
that  this  analysis  is  focused  on  the  movement  of  the  main  surface 
oil  aggregations;  subsurface  plumes  of  oil  may  have  penetrated  to 
peninsular  Florida’s  west  coast  (Paul  et  al.,  2013).  Our  decay  rate  is 
based  principally  on  the  surface  crude  oil  evaporation  rate  (Belore 
et  al.,  201 1 )  -  simulated  subsurface  oil  would  require  a  different 
parameterization. 

The  observed  transport  of  oil  around  the  MRD  into  the  Louisiana 
Bight,  and  then  ultimately  shoreward  to  Grand  Isle  was  well  cap¬ 
tured  by  the  COAMPS-based  oil  forecast.  Coastal  Louisiana’s  com¬ 
parative  misfortune  was  due  not  merely  to  its  proximity  to  the 
DWH  site,  but  also  due  to  the  sustained  surface  conduit  provided 
via  a  buoyancy-driven  current  along  the  Southwest  Pass.  Due  to 
the  unique  cross-shelf  geomorphology  of  the  MRD,  there  is  very  lit¬ 


tle  distance  between  the  MRD  and  the  open  Gulf  of  Mexico.  Indeed, 
the  shelf-break  between  the  MRD  and  the  Mississippi  Canyon  may 
serve  as  an  important  area  of  cross-shelf  water  mass  exchange. 
Both  simulations  and  observations  of  the  DWH  oil  trajectories  sug¬ 
gest  this  is  the  route  the  oil  took  to  transgress  the  outer  continental 
shelf  (50-200  m  isobaths;  see  Fig.  1)  and  precipitate  a  substantial 
landfall  of  oil  along  Louisiana’s  coastline. 

This  finding  may  be  critical  to  understanding  future  distribu¬ 
tions  of  potential  oil  spills  in  the  Gulf.  In  general,  currents  over 
the  continental  shelf  (<200  m  depth)  have  a  tendency  to  flow  along 
isobaths  (alongshore)  and  deep-ocean  properties  are  constrained 
from  transgressing  the  continental  shelf,  as  predicted  by  Taylor- 
Proudman  theory  (see  Brink,  1998,  and  also  Weisberg  and  He, 
2003).  Identification  of  specific  areas  and  mechanisms  that  permit 
‘open  Gulf  to  ‘shelf  water  mass  exchanges  is  required  to  anticipate 
the  fate  of  significant  oil  spills  in  the  major  extraction  region  of  the 
northern  Gulf  of  Mexico.  To  date,  areas  in  the  Gulf  where  the  meso¬ 
scale  circulation  impinges  on  the  shelf  and  the  region  around  the 
MRD  appear  to  constitute  significant  areas  of  open  ocean-to-shelf 
water  mass  exchange  (Biggs  and  Muller-Karger,  1994;  Weisberg 
and  He,  2003;  Jolliff  et  al.,  2008).  Note  that  accurate  forecasting 
would  then  require  regional-scale  knowledge  of  the  circulation 
(Loop  Current  and  associated  eddies)  as  well  as  local  scale  knowl¬ 
edge  of  freshwater  discharges  in  the  MRD  and  potentially  other 
sources. 

These  simulations  did  not  consider  oil  as  a  distinct  surface  layer 
capable  of  responding  to  wind  stress  independently  of  the  ocean 
surface.  Other  simulation  efforts  have  attempted  to  consider  this 
behavior  explicitly  (Sobey  and  Barker,  1997;  Zelenke  et  al.,  2012). 
Nonetheless,  the  qualitative  agreement  between  our  simulated 
17  May  2010  regional  oil  distribution  and  the  17  May  satellite  data 
(Figs.  7a  and  9b)  suggests  that  this  surface  layer  effect  may  be  less 
critical  when  dealing  with  mesoscale  magnitude  oil  spills  in  the 
open  ocean.  It  is  not  well  known  how  the  sea  state  in  the  open 
ocean  will  modify  surface  oil  slick  trajectories  given  the  potential 
mechanical  disruption  of  the  oil  slick,  particularly  at  micron-scale 
thicknesses.  There  is  some  evidence  that  explicit  wind-on-oil 
parameterizations  may  not  be  required  away  from  sheltered  bays 
and  harbors  (Huntley  et  al.,  2011).  Other  studies  seem  to  suggest 
explicit  wind-on-oil  considerations  are  indeed  requisite  (Sobey 
and  Barker,  1997;  Le  Henaff  et  al.,  2012).  A  more  comprehensive 
modeling  treatment  would  require  more  detailed  knowledge  of 
how  oil  of  varying  surface  thicknesses  and  chemical  composition 
will  respond  to  wind  forcing,  sea  state,  and  the  three-dimensional 
hydrodynamic  field. 

Wind  terms  notwithstanding,  both  the  COAMPS-based  and  re¬ 
gional  oil  spill  simulations  presented  here  support  the  assertion 
of  Liu  et  al.  (2011a):  in  the  practice  of  oil  spill  modeling,  ocean  cir¬ 
culation  is  fundamental  to  all  else.  A  key  to  both  sets  of  experi¬ 
ments  is  the  simulation  of  hydrodynamic  conduits  that  may 
expedite  the  transport  of  surface  materials  from  the  accident  site 
to  areas  of  particular  concern.  Additional  considerations  are  then 
contingent  upon  the  scales  of  time  and  space  under  scrutiny.  In 
the  COAMPS  experiments,  simulated  landfall  at  Grand  Isle,  Louisi¬ 
ana  was  accelerated  by  comparatively  swift  coastal  currents 
(~1.3  m  s_1)  contrasted  against  a  background  of  much  more  nom¬ 
inal  surface  current  velocities  (~0.2  m  s_1).  Over  a  distance  on  the 
order  of  ~100  km,  the  timescale  of  transport  for  materials  captured 
by  this  current  is  ~21  h.  On  a  regional  scale,  the  Loop  Current  pre¬ 
sents  a  similar  velocity  hydrodynamic  conduit  for  surface  materials 
(~1.2  m  s_1).  However,  a  consideration  of  Loop  Current  transport 
of  surface  materials  from  the  northern  Gulf  to  the  Florida  Keys 
and  beyond  increases  the  transport  length  scale  by  an  order-of- 
magnitude  (~10001<m).  The  associated  transport  timescale 
(~10days)  is  now  more  commensurate  with  our  estimate  of  the 
half-life  for  surface  crude  oil  (~16  days;  Eq.  (4)).  Thus  weathering 


98 


J.K.  Jolliffet  al./Ocean  Modelling  75  (2014)  84-99 


concerns  become  much  more  pertinent  to  surface  oil  forecasts  with 
the  increase  in  transport  time/space  scales. 

A  remaining  uncertainty  in  this  discussion  is  the  usage  of  dis¬ 
persants.  Over  6  x  106L  of  dispersants  were  released  during  the 
DWH  event  (Judson  et  al.,  2010).  These  materials  are  designed  to 
break  up  the  hydrocarbons  so  as  to  accelerate  weathering,  biodeg¬ 
radation,  and  physical  dispersion.  A  key  remaining  question  is 
whether  or  not  dispersants  were  applied  in  sufficient  quantities 
to  significantly  modify  the  scaling  analysis  presented  above.  We 
note,  however,  that  by  modifying/eliminating  the  buoyancy  restor¬ 
ing  term  and  increasing  (r)  in  Eq.  (5),  our  simulations  may  be  able 
to  provide  an  upper-limit  estimate  of  dispersant  effectiveness  and 
mitigation. 

In  conclusion,  we  have  presented  a  proof-of-concept  oil  spill 
transport  forecasting  method  based  on  the  BioCast  system  and  in¬ 
put  data  from  operational  ocean  circulation  models  and  satellite 
imaging  of  the  ocean.  Given  that  offshore  drilling  will  continue  in 
the  northern  Gulf  of  Mexico  for  the  foreseeable  future,  it  is  proba¬ 
ble  that  oil  spills  of  some  magnitude  may  occur  again.  Shorter  term 
(out  to  96  h)  surface  oil  spill  forecasting  -  with  particular  emphasis 
on  potential  landfall/beaching  of  large  oil  slicks  -  is  critically 
dependent  on  accurate  ocean  current  forecasts  and  knowledge  of 
where  cross-shelf  water  mass  exchanges  are  likely  to  occur.  In 
our  particular  example,  this  cross-shelf  exchange  is  critically 
dependent  on  accurate  shoreward  fluxes  of  buoyancy.  Longer-term 
simulations  for  oil  slick  transport,  likely  required  when  oil  spills 
are  of  regional  scale,  need  to  more  carefully  consider  the  intrinsic 
dynamics  of  oil  weathering  processes  and  potential  oil  source 
terms.  By  considering  both  the  timescales  of  the  degradation  pro¬ 
cesses  in  concert  with  material  transport  pathways  driven  by  the 
ocean  circulation,  our  simulations  did  not  indicate  any  significant 
surface  oil  contamination  beyond  the  northern  Gulf  of  Mexico. 
Computer  simulations  used  in  the  future  for  oil  spill  response  must 
consider  the  timescales  of  all  the  processes  involved. 
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